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temperature fluctuations in the source of CR with the "knee" structure reflecting an abrupt change of the 
pattern of such fluctuations. This would result in a generalized nonextensive statistical model for the 
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the resulting chemical composition of the CR, which follows the experimentally observed abundance of 
nuclei. 
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1. Introduction 

The energy spectrum of cosmic rays (CR) has characteristic power-like behavior with a "knee" structure (plus 
some other less prominent features) and remains constantly matter of hot debate (see [1] and references therein) . It 
could reflect the action of different regimes of diffusive propagation of CR in the Galaxy combined with its different 
chemical composition, but it could also be due to some, so far unspecifled, property of the production processes 
within the source of the CR itself. In this work we shall consider this possibility assuming that CR are produced 
following a generalized nonextensive thermal approach [2, 3]. Actually, nonextensive statistical mechanics [2] has 
been applied to CR before: in [4] the "knee" structure was attributed to the crossover between the assumed fractal- 
like thermal regimes of CR propagation (characterized by different temperatures T and nonextensive parameters 
q), whereas in [5] the possible nonextensive thermal features of CR flux have been investigated but only up to 
the "knee" region, the origin of which was not discussed. In both cases the obtained values of temperatures were 
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much too high to be accommodated by any known physical mechanism. In this paper we propose a mechanism 
which is apparently capable to describe the whole spectrum of CR, including the "knee" region, using physically 
reasonable values of temperature of the source of CR. The observed power-like behavior of the energy spectrum 
of CR is attributed (as in [5]) to fluctuations of the temperature in the source producing CR and the occurrence 
of "knee" (cf., Fig. la) is connected with some abrupt change of this fluctuation pattern [6], visualized by a 
dramatic change in the nonextensivity parameter observed Fig. lb. However, to keep the temperature of the CR 
source acceptably low, one has to allow additionally for energy transfer to the production region; this is assumed 
to proceed through the mechanism proposed in [3] and is characterized by some effective temperature Te// ^• 
These points summarize what we call a generalized nonextensive thermal approach (GNTA), which we shall now 
describe in more detail. 

It must be stressed at this point that, for the sake of clarity of presentation, we consider in what follows only a 
very simplified situation. Namely, we assume that GNTA is, for a moment, the only mechanism of production 
of the CR present. It must be realized that in reality GNTA would have to be incorporated into many other 
possibilities considered in the usual analysis of CR (and listed, for example, in [1]). 

The organization of our paper is as follows. In the next Section we present some basic considerations concerning 
nonextensive statistics and CR, out of which the discussion of the chemical composition of CR seen from that 
point of view is a new element here. Section 3 contains our results and their physical interpretation in terms 
of some specific properties in the superfluid stages of neutron stars supplied by the proposition of introducing 
phenomenologically the energy transfer to CR (described by some effective temperature T^ff and needed to assure 
the consistency of obtained parameters). In Section 4 we discuss the influence of acceleration and propagation of 
CR on energy spectra and composition. A summary and concluding remarks are presented in Section 5. 

2. Basic elements of nonextensive statistics and Cosmic Rays 
2.1. Generalities 

Nonextensive statistical mechanics as proposed and developed in [2] is based on the generalized entropy functional 
(Tsallis entropy), 

_ SdEP^jE) - 1 
^« - —I ■ 

Its maximization under appropriate constrains yields a characteristic power-like distribution (g-exponential dis- 
tribution, expg(. . . )): 




(2) 



^ Actually, this mechanism was originally invented to describe some features of heavy ion collisions [7] in which 
energy was transferred out of the system; tn the case of the CR it is transferred towards the system. 
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For q 1 one recovers the usual Boltzmann-Gibbs-Shannon (BGS) entropy and the usual exponential distribution. 
This equilibrium distribution can alternatively be obtained by solving the following differential equation, 

dP{E) _ P^jE) 
dE ' T ■ ^■^> 

The extended version of this equation with two terms (accommodating two different values of q and T) has been 
used with apparent success in [4] to describe the flux of CR. The "knee" appears there as a crossover between two 
fractal-like thermal regimes characterized by (T, g) and {T',q'). However, the values of temperatures obtained 
there {{T,T') ~ 100 -i- 1000 MeV) are uncomfortably high and cannot be attributed to any known mechanism of 
CR production. 

On the other hand, there is growing evidence that a nonextensive formalism applies most often to nonequilibrium 
systems with a stationary state that possesses strong fluctuations of the inverse temperature parameter (3 = 1/T 
[6, 8] (cf., also [9]). In fact, fluctuating /3 according to gamma distribution (cf. Eq. (9) below) with variance 
Var{/3) results in a power like distribution (2) with the deviation of the nonextensivity parameter q from unity 
being given by the strength of these fluctuations, 

This observation was used in [5] to describe the energy spectrum (but only up to the "knee" region). Again, 
although the results were reasonably good the estimated temperature T ~ 170 MeV is far too high. This was 
because the author insisted on the description of the whole range of energy spectrum up to the "knee" region, 
including its very low energy part, which is, however, governed mainly by the geomagnetic cut-off and diffusion 
effects and should therefore be considered separately (it is thus not covered by our approach). 

2.2. Energy spectrum 

Treating CR as relativistic particles (for which rest mass m can be neglected) their energy is _B ~ p and the 
density of states is that of an ideal gas in three dimensions, Q{E) oc E^ . The corresponding energy spectrum 
$(£) is then 

^{E) = NoE^P{E), (5) 
where A'^o is normalization factor. For P{E) given by Eq. (2) we have, for E » T, power spectrum 

'f{E)o,E-^; 7=^=^. (6) 

As seen in Fig. la it changes in the region named "knee" where the slope 71 ~ 2.7 at energies below ~ 10^^ eV 
and 72 — 3.1 above it. In the language of the nonextensivity parameters it would mean that qi = 1.213 before 
and §2 = 1.196 after the "knee". For q understood as a measure of fluctuations, as it is the case in our paper, one 
therefore witnesses at the "knee" a change of fluctuation pattern. 
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Fi gUrG 1. (Color online) (a) Schematic view of the observed CR energy spectrum [1] with characteristic features clearly indicated. 

(b) The pattern of fluctuations obtained from the above by means of nonextcnsivc statistics approach and represented by 
the heat capacity Cv as given by Eq. (11) (see text for details). Notice that previous irregularities arc now dramatically 
enhanced. 



2.3. The chemical composition 

If the energy distribution of CR follows the Tsallis formula (2) , it is natural to expect that also the abundance of 
nuclei with mass A will follow the same pattern. We therefore expect that 

Lj{A) oc 

where e is the average energy of nucleon equal to e = |-Ef. The typical value of the Fermi energy for nucleus 
consisting of A nucleons (distributed in sphere of radius R = 1.25^1^''^) is Ef — 30 MeV, it means then that 
e ~ 22 MeV. In Fig. 2 we show u;{A)/ij{A = 1) for e = 22 MeV and T = 100 MeV as function of A for two values 
of g: q — 1.2 and q = 1.15. As one can see the sensitivity to q is rather weak ^. The average mass number (for 
the spectrum i^{A)) is 

Numerically evaluated {\n{A)) equals 1.82 below the "knee" (for greater q) and 1.78 above the "knee" (for smaller 
q) and shows that predicted changes of chemical composition (due to changes of spectral index) in the "knee" 
region are negligible. Notice that the upper limit estimation for the usual BG distribution (corresponding to 
q = 1 here) leaves the majority of points for large values of A well above the curve. From this point of view our 
prediction is much better (although, the thermal model alone is already able to provide quite reasonable chemical 

^ In fact, the relative abundance shown in Fig. 2 depends on the temperature chosen. However, reasonable fits are 
possible only in the limited temperature interval, 60 Me V < T < 160 Me V, and for the nonextensivity parameter 
satisfying roughly the relation q — 1.5 — 0.29r/100. Our fit ts for the same parameters {T,q) as used to describe 
the energy spectrum . 



1-(1-, 



As 



(7) 
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composition of CR) 

One should keep in mind that the relative dissemination of nuclides shows many intriguing features which should 
be connected with specific properties of nuclei, cf., for example, [11]. We close this section by stressing that our 
io{A) as given by Eq. (7) (and therefore represented by a continuous curve) does not account for any differences 
between nuclei. In fact, it even does not describe the abundance in the source. It describes only the emissive 
power of thermal source and shows that this factor (neglecting differences between particular nuclei) determines 
the global characteristics of the observed spread of nuclides. 

3. Results 

3.1. Temperature fluctuations 

As mentioned above, the special role in converting an exponential distribution to its g-exponential counterpart 
play fluctuations of the inverse temperature /3 described by a gamma function [6, 9], 

= ^(M/3)"'-'exp(-/3/.) (9) 

A remark of caution is in order here. Our {A) agrees with observations in the low energy region but the observed 
composition shows changes with energy. On the other hand, we show small changes connected with the change of 
the spectral index 7 and this could be caused by some other mechanism influencing chemical composition which we 
have not accounted for. The most important is problem of energy dependence of the effective temperature T^f f ( cf. , 
Eq. (18)). In fact, we do not know To and cannot say (numerically) how Tvisc changes with the energy. Instead, 
we just put roughly Tej / = 100 Me V. A more exact analysis of experimental data, including energy dependence of 
the chemical composition, would be helpful in estimation of To itself. We plan to address this problem elsewhere. 
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where jj."^ = [3o{q — 1) and = q — 1. There are a priori at least two scenarios leading to such /(/3): (i) one 
can have many sources with different temperatures, the number of which is distributed that way or (ii) one can 
have temperature fluctuations in small parts of a source. The first possibility is, however, rather unlikely because 
in this case either one would have to accept sources with unphysically large and small temperatures or else use 
the temperature distribution in some limited domain, i.e., work with a truncated version of gamma distribution. 
This, however, would result in a very characteristic rapid break in the energy spectrum. This is not observed in 
the experiment. 

We shall therefore concentrate on the second possibility. To illustrate it, suppose one has a thermodynamic 
system, different small parts of which have locally different temperatures, i.e., its temperature understood in the 
usual way fluctuates. Let £^{t) describes the stochastic changes of temperature in time and let it be defined by the 
white Gaussian noise ({^(t)) ~ and {(,{t)^{t + At)) = 2DS{At)). The inevitable exchange of heat which takes 
place between the selected regions of our system leads ultimately to an equilibration of temperature and, as shown 
in [6], the corresponding process of heat conductance eventually leads to the gamma distribution (9) mentioned 
before with variance (4) related to the heat capacity Cv of this system (expressed in units of Boltzmann constant 
fcs, which we put equal unity in what follows) by [12] 



and we have finally that 



Cv = ^=7 + 2, (11) 
q-1 



where we have used Eq. (6) connecting the spectral index 7 of the energy spectrum with the nonextensivity 
parameter q. In this way, we come directly to the possible physical interpretation of the nonextensivity parameter 
which allows us to translate the pattern observed in Fig. la into that shown in Fig. lb. Here Cv, as given by 
Eq. (11), is shown as a function of energy at which we observe the essential influence of given Cv on the slope 
of energy spectrum. As will be discussed in detail below in Sec. 3.2, although Cv as such is energy independent 
it can change with temperature T in the source and changes abruptly at some temperature Tc (identified, for 
example, with Tcut in Sec. 3.2). These changes of Cv result in changes of the slope of the observed energy 
spectrum at i5 ~ 2Tc). 

To summarize this part: In what follows, we shall concentrate mainly on the change of fiuctuation pattern in 
the "knee" region and we shall argue that it could indicate an abrupt change in the heat capacity in the source 
of CR of the order of C2/C1 — 1.09. Notice that this change is much more pronounced and dramatic than the 
corresponding change of slope in the "knee" region observed in Fig. la. 

3.2. Possible physical interpretations of the fluctuation pattern 

Can one expect something of this kind to happen in the astrophysical environment of the CR? In what follows we 
shall argue that, indeed, one can. Let us first notice that the subject of temperature fluctuations in astrophysics 
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is a much-discussed problem nowadays. Its eflect on the temperatures empirically derived from spectroscopic 
observations was first investigated in [13] whereas in [14-17] it was shown that temperature fluctuations in 
photoionized nebulea have great importance to all abundance determinations in such objects. This means that 
discussion of the heat capacity or, equivalently, the behavior of the parameter q defining the energy spectrum, is 
fully justified. In what follows we shall concentrate on the problem of the heat capacity of astrophysical objects, 
in particular in neutron stars, concentrating on some peculiarities connected with their description in terms of 
Fermi liquids. 

We start with fluid/superfluid transitions in such systems and their effect on the heat capacity. In neutron stars 
one observes the following feature. The total specific heat of their crust, C, is the sum of contributions from the 
relativistic degenerate electrons, from the ions and from degenerate neutrons. In the temperature that can be 
reached in the crust of an accreting neutron star (which is of the order of T ^ 5 ■ 10* K and is below the Debaye 
temperature Td ~ 5 • 10^ K) we have don < Ce < Cn- When the temperature drops below the critical value 
T — Tc the neutrons become superfluid and their heat capacity C!^^ increases [18, 19], 



— ^3.15-exp 



-1.76 



T 



2.5-1.66 I I +3-1 



T 



(12) 



At T ~ 0.7Tc we have C^^ ~ 1.1C„ what corresponds to the changes of spectral index by A7 ~ 0.5. To 
summarize: one witnesses here an abrupt change in the heat capacity at some temperature, i.e., a phenomenon 
we were looking for. 

The above example tells us that it is reasonable to expect the abrupt change of the specific heat in the CR source. 
Assuming that this really happens and taking seriously the apparent connection between Cv and q expressed 
by Eq. (11), we are lead to the natural conjecture that in this case the usual fluctuation pattern given by the 
gamma distribution (9) should be modifying accordingly. It is then assumed to be done by replacing Eq. (9) by 
its slightly modified version, characterized by two nonextensivity parameters, qi (acting before some temperature 
Tcut) and 52 (acting after Tcut)- The change qi — » g2 at Tcut is assumed to be abrupt and the temperature Tcut 
becomes a new parameter in our description. Following our proposition one obtains the following flux of CR: 



${E) = iVo-B 



l/Tc, 



No ■ [P,, (E) ~ a,{E)P,, (E) + a2{E)P,, (E)] , 



(13) 



where Pq. (E) are given by Eq. (2) and 



1 1 - (1 - gQ E/T 
gi — 1 ' {qi — 1) Tcut/T 



(14) 



Our results are presented in Fig. 3 where q = 1.214 in Fig. 3a whereas in Fig. 3b gi — 1.214 and g2 = 1.2; in both 
cases T — 100 MeV. Notice that now we do not have a spectrum where, as in Section 2.2, we change the value of 
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Figure 3. (Color online) (a) CR energy spectra fitted by the single Tsallis distribution with q = 1.214 and T = 100 MeV. 

Notice the total inadequacy of a simple exponential (Boltzman, denoted by BG here) distribution in description of 
these spectra, (b) The CR energy spectrum fitted by the double Tsallis distributions discussed in text, cf. Eqs. (13) 
and (14). It was obtained for T — 100 MeV and assuming that fluctuations of the temperature change abruptly at 
Tcut Ri 10^^ eV Ki 10^" K from q = 1.214 to g = 1.2. Data are from [20], cf. also [1]. 



q at some energy Ec to get the observed structure. Spectrum (13) is obtained by changing q at some temperature 
Tcut (i.e., by changing slightly the shape of gamma function (9)), this means that each q gives a spectrum for 
all energies. For this reason the parameters q here have slightly different values from these in Section 2.2. With 
spectrum given by Eq. (13) the "knee" region is reproduced very well, however, the price to be paid is the need 
of a suitable choice of temperature at which the fluctuation pattern changes (which amounts to assume the value 
of Tcut « lO''^ eV ^ lO^'-* K) 

The above mechanism is only able to describe the "knee" region. To describe all details seen in Fig. lb let us 
consider an other feature of heat capacity in Fermi liquids, namely its dependence on the effective mass of nucleons 
consisting such liquid. Following [18] the proton heat capacity is proportional to the ratio of the effective mass 
of the proton in the neutron fluid to the mass of the free proton, C ^ m* /m. In the case of a mixture of Fermi 
liquids the proton effective mass m* is affected by interactions with neutrons and other protons and is given by 



m 3 



fVV I / \ rpn 



(15) 



where Dp denotes the density of quasiparticle states at the Fermi surface given by wave vectors kp^ and kpp 
for, respectively, neutrons and protons, whereas /f and ff" are Landau parameters [26]. Fig. lb can then be 



Two remarks: (i) The nucleon superfluidity was predicted already in [21] and today pulsar glitches provide strong 
observational support for this hypothesis [22]. Nucleon superfluidity arises from the formation of Cooper pairs 
od fermions (actually m [23] also quark superfluidity from cooling neutron stars were investigated). Continuous 
formation and breaking of the Cooper pairs takes place slightly below T = Tc ( critical temperature Tc is in the 
order 10® — 10^" K). (ii) Neutron stars are born extremely hot m supernova explosions, with interior temperatures 
around T ~ lO^'^ K. Already within a day, the temperature in the cental region of the neutron star will drop down 
to oc lO'' — 10^*^ K and will reach 10^ K m about 100 years [24-]. The first measurements of the temperature of 
a neutron star interior ( core temperature of the Vela pulsar is T ^ 10* K, while the core temperature of PSR 
B0659+14 and Cemmga exceeds 2 • 10** K) allow us to determine the critical temperature Tc ~ 7.5 • 10® K [25[. 
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interpreted as showing changes of C with energy in the Fermi hquid. We start with the superfluid hquid with 
Ci =4.7 (here m* represents effective mass for pp and pn interactions), when energy increases we stop to see 
nuclear interactions and C2 = 5.1 (with m* representing pp interactions only), finally, for large T, one has the 
Fermi gas with C3 = 5.3 and, still further, the usual Fermi liquid ^ . Notice that 

this results in the following relation between Landau parameters. 



fr C2 - C3 



(17) 



In the case of a one-component Fermi liquid we have the well known identity, m* /m = 1 + Ff'^ /3, where 
= Dp/f . From (17) we see that in a two-component Fermi liquid the quantity 1 — m*/m is 3 times bigger 
(this is because the parameter /i which determines the interaction between quasiparticles is negative, resulting 
in smaller effective mass). From properties of excited states in nuclear matter {Pb and neighboring nuclei [28]) 
F^" = -0.5 ± 0.25. If F^ < F"" < Ff^ and taking (after [29]) Ff" - F^^ = -0.2, we can estimate that for 
neutron-star matter one has m* /m = 1 + ~ 1 — 0.4 ± 0.3 = 0.6 ± 0.3. 

3.3. The notion of the effective temperature T^ff 

Let us now come back to the results presented in Fig. 3. Although a double Tsallis fit looks rather impressive 
there are two shortcomings which we shall now discuss in more detail. First is the fact that we still need a too high 
value of the temperature, T = 100 MeV, which cannot be accommodated by any reasonable physical mechanism 
of production of CR's. Second is the very high value of the T^ut temperature where change in the fiuctuation 
pattern is supposed to take place. 

The possible way out of both dilemmas we are going to propose now is to keep the value of T = 100 MeV but 
change its meaning. This can be done by adding to the mechanism proposed in [6] (which was accounting only for 
the possible fluctuations of T) an additional effect of the possible viscosity which describes the possible transfer 
of energy between the region of production and surroundings (cf. [3] and [7] for details). As a result one gets the 
same power-like distribution as before but with the previous T replaced by an effective temperature T^ff : 

T,ff =To + iq-l)n^sc■ (18) 



^ It is worth to remember that fluctuations of temperature we are talking about in this work refer to fluctuations 
in a small region V . For a Fermi liquid the heat capacity expressed tn umts of Boltzmann constant ks (i.e., for 
kB ~ 1) is of the order C ~ 3 • 10"^^ cm~^ [27]. Therefore, taking values of C estimated from the slope of the 
primary CR spectra (cf. Fig. 3) one gets that the size of the region of fluctuations is V ~ 10* frr? . 
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Here To is the temperature around which one has fluctuations and T^isc is some new parameter depending on the 
transport properties of the surrounding space around the emission region. 

A few words of explanation are necessary at this point. Following [7], in the case of heavy ion collisions, where 
this concept has been introduced for the first time, T^isc = rif{u)/{Dcpp), where rjf{u) presents the effect of a 
possible viscosity, with viscosity coefficient -q and f{u) = (^fj:^ + being velocity), whereas D, Cp and 

p are, respectively, the strength of the temperature fiuctuations, the specific heat under constant pressure and 
the density. Here this quantity is supposed to model the possible transfer of energy towards the CR particle. To 
estimate its value let us notice that (g — l)Tvisc ~ {t/tc)To where r and Tc are, respectively, a relaxation time of 
the corresponding dissipative process and the mean collision time. They have distinct physical meanings. Namely, 
T is the time (in most cases macroscopic) taken by the system to spontaneously return to the steady state (whether 
in thermal equilibrium or not) after being suddenly removed from it. It is, to some extent connected to the mean 
collision time Tc of the particles responsible for the dissipative process. There is no general formula linking them, 
their relationship depends in each case on the system and circumstances under consideration. In particular, r is 
small for photon-electron and photon-photon interactions at room temperature (oc 10^^^ and oc lO"'^^ seconds, 
respectively). In the case of neutron star, r is of the order of the scattering time between electrons (which carry 
the heat) but this fact is not obstacle (no matter how large the mean free path of this electron may be) to consider 
the neutron star as formed by a Fermi fluid of degenerate neutrons. One finds in this case the relaxation times as 
large as r ~ 10^ seconds (for T ~ 10^ K, p ~ lO^^g/cm'^ and u ~ 10^ cm/s) [30] whereas in the degenerate core of 
aged stars r can reach 1 second [31]. The fact that r can quantitatively greatly differ from Tc is most dramatically 
suggested by the matter-radiation decoupling in the early universe (the relaxation time of shear viscosity turns 
out to be several order of magnitude larger than the collision time for most of the radiative era) [32]. For this 
region one can expect that Tvisc may be large enough (for example, Tvisc/Tq ~ 10*) to be a reasonably parameter 
describing the influence of surroundings on the CR particle. Therefore in such picture there are fluctuations of 
temperature T characterized by the parameter q and connected with the character of the source, for example, 
as discussed here, but there is also a surrounding space around the emission point which can pump energy into 
some selected region from which a CR particle is emitted. This is described by the positive parameter T^isc- As 
a result we are getting the total Te// that grows with q. 

Our schematic view in what concerns the fate of the CR spectrum is then the following, (i) CR are produced 
in an object where we have low temperature which experiences some fluctuations around T = To given by g > 1 
(with To of the order of MeV). Now To can be as small as desired to properly fit the expected features of the CR 
source, (ii) By introducing to the previous description [6] some (so far unspecified) flow [3, 7], the corresponding 
Tsallis distribution now has T — Tefj, which can be quite large for example, of the order of the observed 100 
MeV). (iii) Therefore, in the spectrum of CR one observes a Tsallis distribution with Tejf ~ 100 MeV (as in the 
previous attempts [4, 5]) but now it is not the temperature of the source of CR itself, but it is composed from 
the temperature of the CR's source itself. To, and the effect of the action of the surrounding space which is given 
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by Tvisc- It should be stressed that such an effect appears only in the case of Tsahis distribution. Only for g > 1 
is Te// > To- If there are no temperature fluctuations one always has T^ff — Tq. 

Interestingly enough, the proposed mechanism seems to be capable also to alleviate the second problem, namely 
the uncomfortably high value of the Tcut where the change of the nonextensivity parameter q occurs. This is 
because, as discussed in [43], when one considers viscosity effects in stars, one observes very large abrupt changes 
of the viscosity coefficient which can easily result in an effective T^isc of the desired order to fit the parameter 



4. Propagation and sources of cosmic rays 

Following [4] one can argue that nonextensive approach represented by Eq. (2) can also be connected with the 
process of CR propagation. In our case such scenario would mean that parameter Tvisc in Te// can represent 
summarily some effective modelling of the propagation of CR. Whereas detailed discussion of such a possibility 
is outside the scope of the present paper, a few words on propagation of CR and its possible influence on the 
spectra and composition of CR are in order here. 

So far we have considered only the sources of CR insisting on the possibility that they are of thermal but nonex- 
tensive type. We shall now address the possible influence on our results by the propagation of CR. Measurements 
of the composition and energy spectra of CR characterize the CR population after the observed particles have 
travelled from their production sites (i.e., the sources) through distant space towards the remote detectors. It is 
expected that both the composition of the particles and the shape of their energy spectra undergo in such process 
changes due to a variety of processes encountered during the propagation. 

The results presented in this work concern only production in the source, which is supposed to be thermal and 
nonextensive for some reasons discussed here. However, the fact that one can confront them with experimental 
data indicates the possibility that they depend only weakly on the propagation process in what concerns the 
slopes of the energetic spectrum and the global characteristics of the chemical composition (notwithstanding the 
fact that propagation is very important factor). 

The possible way out of this dilemma is to argue that the break in the original spectrum and connected with it 
phase transition occurs actually at much slower energies and that resultant spectrum is then accelerated to the 
observed energies - for example by the magneto-hydrodynamical turbulence and/or shock discontinuities (i.e., by 
the so called diffusive shock acceleration (DSA) mechanism, cf. [33]). The simplest version of this mechanism, as 
discussed in [34], implies that distribution function after the shock, fafter{p), is related to the original distribution 
before the shock, fbeforc{p) in the following way: 



where b = 3r/{r — 1). Here p denotes the particle momentum, r = P2/P1 describes the compression of densities 
across the shock and pmin denotes the minimal value of momenta. DSA mechanism transforms a 5{p — po) 



Tcut • 




(19) 
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spectrum of relativistic particles in a power-like spectrum of the type f{p) oc p^*". For example, if one has initial 
spectrum of the form oc which encounters a shock with strength given by b then Eq. (19) shows that: (i) in 
the case of 6 > c (which corresponds to the initial spectrum injected into the shock being softer than it would 
result from a (5-function) one has faftcr oc p^'^ (i.e., the acceleration process does not change the shape of the 
spectrum); {ii) in the case of 6 < c (i.e., for the steep initial spectrum) one has f after oc p~^ , which coincides with 
result of injection of a (5-function into the shock. It can be then shown from Eq. (19) that when the strength 
of shocks is larger then the slope of the injected spectrum, the shape of the spectrum should be given by the 
production spectrum in source [34]. 

To summarize this part: stochastic mechanisms of acceleration of CR particles (like acceleration on the fronts of 
shock waves or Fermi acceleration in turbulent plasmas, both analogous in some sense to Brownian motion) do not 
change the shape of the power-like production spectra, but, unfortunately, they are not particularly effective, i.e., 
they do not lead to large increase of energy [35]. The increase of energy per one collision is of the order / E ~ 
/c^ , what for the plasma velocity u ~ lO" — 10^ cm/s gives AE/E ~ 10~* and leads to the mean relative increase 
of energy during the time life of Galaxy {t ~ 3-10^^ s) only by factor {S) ~ 3. Fluctuation on the steep spectrum of 
accelerated particles result in additional increase of energy. Because of the multiplicative character of acceleration 
we have log-normal distribution of variable S, P{S)dlnS = (\/27r(7) exp [- (In 5 - ln(<5))^ /2(t2] din (5, what 
results in shift of the spectrum of source on the energy scale by {S)^~'~^^^~' exp [(7 — 1)^(7^/(27)], where is 
variation of the distribution P{5). For j ^ {S) ^ a'^ ^ 3 we can obtain only order of magnitude shift of the energy 
spectrum. 

Diffusive propagation of CR component is commonly summarized in a continuity equation for the differential 
density, Ni{E), of each component [36]. If convection effects are neglected (which is probably valid approximation 
at high energies) together with effects due to energy gain or loss and to radioactive decay, the continuity equation 
becomes [37] 

1 [DJPA AT. 1 

(20) 



1 

A-^E)+As'iA) 



3,(S) iV, 

— h 2jfc>l -7 

pep Ak^i 



Here Qi{E) is the rate of production in the source and A^^i quantifies the probability of a nucleus k to spallate 
into a product i in an interstellar interaction. The two quantities, the propagation length, A{E) — /3cpT{E), and 
the average spallation length, As{A) = m/a{A), characterize the propagation of cosmic rays and the change of 
the atomic number (the loos of nuclei) due to the spallation effect (here f3 is velocity, p the density of the material 
in the galactic space and t{E) is the average time which particle spend in the galaxy). 

The propagation path length decreases with energy but it is assumed that it has the same value for different 
nuclei of the same rigidity. On the other hand, the spallation path length depends on the atomic number A 
(essentially like oc A^^''^), its energy dependence remains, however, very weak for the relativistic particles 
and is therefore neglected. For high energies (above TeV/amu) A < As and for the approximate scaling behavior 
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Ak^i = Ak-*i (Ai/Ak) the relative abundance is roughly given by 

jV(A,g) QiA,E) 

N{Ak,E) Q{Ak,E)- ^ > 

For stable nuclei the observed abundances do not difTer substantially from the relative abundances in the source 
[37]. 

Some remarks concerning the applicability of Eq. (20) are in order here. It is frequently assumed that the 
propagation path length A decreases as function of energy, A oc E^"''' [38]. Since the interaction length is almost 
independent of the primary energy this necessitates Q{E) oc E~^'^ spectrum at the sources to explain the observed 
N{E) oc E^"^'"^ spectrum at the Earth. Because of this, the values of q evaluated by us from fits to the observed 
spectrum must be regarded as the corresponding lower limits of the nonextensivity parameter q {q m 1.213 for 
7 = 2.7 and q ~ 1.244 for 7 = 2.1). In fact, because of large uncertainty in what concerns the energy dependence of 
the propagation path length which still exists, we cannot present exact values of q in the source. For example, the 
model presented in [39] predicts flatter than the above mentioned spectra at the sources before the knee requiring 
therefore a stronger dependence of A{E) on energy. Recent measurements of the TeV gamma ray flux from a shell 
type supernova remnant yield spectral index 7 — 2.19 ± 0.09 ± 0.15 [40], in agreement with the standard model 
^. On the other hand, for the Crab Nebula a steeper spectrum (with 7 = 2.57 ± 0.05) has been obtained [42], 
indicating that probably not all sources exhibit the same behavior. Moreover, the A oc E^^''^ dependence of the 
propagation path length cannot be extrapolated to the knee energies whereas taking A oc E~^''^ dependence, as 
discussed in [38], necessitates additional assumptions concerning the spectral shape at the source, Q{E), in order 
to explain the observed spectra with spectral indices in the range 7 ~ 2.55 ± 2.75. 



5. Summary and conclusions 

We have proposed and discussed the possibility that CR can originate from nonextensive thermal sources described 
by a nonextensive formalism proposed in [2]. Our motivation was the observation that the spectrum of CR has, 
in general, a power-like shape, E~^ , and such behavior is naturally accounted for in a nonextensive approach 
with nonextensivity parameter q= (3 + 7)7(2 + 7). Looking more closely one encounters a characteristic "knee" 
structure in this power-like behavior, with 71 ~ 2.7 at energies below Eknee ~ 10^^ eV and 72 — 3.1 above it. 
This can be also accounted for in a nonextensive approach with two different values of nonextensivity parameter: 
qi — 1.213 before the "knee" and q2 — 1.196 above it. 

^ The Standard Model for Galactic Cosmic Rays is based on Supernova Remnant (SNR) paradigm and includes 
four basic elements: (i) SNRs as the sources, {ii) SNR shock acceleration, (Hi) rigidity dependent injection as 
mechanism providing the observed CR mass composition and (iv) diffusive propagation of CR m the galactic 
magnetic fields [4-1 j- 
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From our previous experience with applications of the nonextensivity to different physical processes (cf., [3] and 
references therein) we can trace the origin of such power-like behavior back to some intrinsic, nonstatistical 
fluctuations of temperature in the CR's source [3, 6, 9]. In this case the nonextensivity parameter q is regarded 
as a measure of the heat capacity (see Eq. (11)). This means that the measured energy spectrum (Fig. la) can 
be converted to the energy dependence of the heat capacity C. The result is shown in Fig. lb. As one can see, 
Cv acts here as a kind of magnifying glass converting all subtle structures of ^{E) into much more pronounced 
and structured bumps. Its importance would parallel the long-standing discussion of the origin of the "knee" -like 
structure of the energy spectrum, but exposed in a much more dramatic and visible way. 

As a plausible physical mechanism leading to changes in Cv of the order of C2/C1 = 1.09 (corresponding to 
change in q, describing the observed change in the spectral index 7) we have proposed fluid/superfluid transitions 
in Fermi liquids used to describe neutron stars, which we model by suitably modifying the gamma distribution 
(9) describing temperature fluctuations (essentially by changing parameter q at some temperature Tcut)- To 
get fits as presented in Fig. 3, while at the same time keeping temperature in the CR's source. To, acceptably 
low (of the order of MeV, i.e. of the order of the interior stars temperature) we have to resort to an approach 
allowing not only for fluctuations of temperature but also for the energy transfer to the production region from its 
surroundings introduced recently in [3, 7]. This allows us to keep the critical temperature (corresponding to the 
nucleon superfluidity) around Tc ~ 0.1 — 1 MeV with effective temperature used in the fits remaining as high as 
Tef f = 100 MeV. This also allows to quantitatively understanding that the origin of changes of the nonextensivity 
parameters at the temperature as high as Tcut — 10^^ eV ~ 10^® K required in our approach could be in some 
specific viscous effects in stars [43]. 

One should stress at this point that the mechanism we proposed, namely that CR can indeed originate from 
nonextensive thermal sources, must be, for a while, regarded only as a plausible scheme, which would have to be 
checked together with other mechanisms aiming to describe the CR spectra, their composition and propagation 
[1]. The need for such an analysis (which, however, goes outside the limited scope of this paper) is, for example, 
visible when one realizes the following. On one hand one observes that the energy spectrum of CR depends mainly 
on the nonextensivity parameter q. Dependence on T is visible only for low energies. For E >> T we observe 
scale-free, mostly T-independent behavior. The temperature has therefore marginal influence on the shape of 
the energy spectrum of CR. On the other hand, a quite opposite situation is encountered when considering the 
chemical composition of CR where (A) depends only very weakly on q but it depends linearly on T = Te//, see 
Eq. (8). It is then plausible that analyzing simultaneously the energy spectrum and composition one could obtain 
both q and T^ff (i.e., according to Eq. (18), the To of the source and T^isc responsible for the energy transfer). 
Notice that around the "knee" (where one expects changes in the chemical composition) one can essentially freely 
vary Te// without affecting the shape of the energy spectrum but substantially changing the chemical composition 
(A). 
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